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Abstract 

In this review, we discuss some interesting issues in charm physics which is fuh with puzzles 
and chaUenges. So far in the field there exist many problems which have not obtained satisfactory 
answers yet and more unexpected phenomena have been observed at the present facilities of high 
energy physics. Charm physics may become an ideal place for searching new resonances and 
studying non-perturbative QCD effects, moreover probably is an area to explore new physics beyond 
the Standard Model. More data will be available at BESIII, B-factories, LHC and even future ILC 
which may open a wide window to a better understanding of the nature. 
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I. INTRODUCTION 



The charm quark was a long-expected member of the quark family. It is noted that if 
only u quark is the intermediate fermion at the s-channel, the cross section of the scattering 
s + — > (i + would increase with the incoming energy and is unacceptable in physics. 
It demands existence of another species of quarks having the same charge as u quark which 
serves as an additional intermediate fermion to compensate the bad high energy behavior 
for the process, i.e. retain the unitarity. th,s new speces is the char,n quark fl. Moreover, 
without the charm quark the anomaly in the electro- weak model cannot be cancelled 
so that the renormalizability of the whole theory would be spoiled. Later by studying the 
— mixing, Gaillard and Lee [sl estimated that the mass of the charm quark should 
be around 1.5 GeV. Thus all the urgency of saving the beautiful theory appealed to discover 
this charming "charm" quark. Then the discovery of J/ip meson and other members of the 
if) family became a milestone of particle physics The ground state of the family are J /ifj 
and rjc whose mass is about 3.1 GeV and 2.98 GeV, so roughly it implies that nic ~ 1.5 
GeV which is amazingly consistent with Gaillard and Lee's estimate. But definitely it is not 
the end of story for the quark family, the discovery of bottom quark requires existence of 
its partner which evaded observation for a very long while until it was eventually found at 
the TEVATRON and possesses an astonishing heavy mass about 176 GeV. Till then the 
three generation structure of quarks and leptons seems complete, even though the fourth 
generation of quarks and leptons is still under discussions. 

We know, the m, d and s quarks reside in a triplet of the global SU(3) quark model which 
successfully describes relevant phenomenology, but there is, so far, not a special symmetry 
to associate the rest three heavy quarks^. The charm quark may be a special one in the 
quark family, because it is heavier than the first three light quarks and does not belong to 
the regular flavor SU(3), but stands in a weak doublet with the light strange quark. The 
charm is not light at all, but also not too heavy as the bottom and even top quarks. The 
intermediate mass determines the special characteristics of hadrons which contain charm 



^ There indeed is the so-called Heavy Quark symmetry SUj{2) ® SUs{2) (we will discuss later) to connect 
the b and c quarks, but that symmetry is a symmetry when the quark mass approaches infinity and mainly 
can simplify the calculation of transition from b to c. It is not like the SU(3) for light flavors and is also 
not a symmetry in the common sense. 
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and/or anti-charm. Recently, some researches count that the first four quarks are different 
from the last two heavy ones: bottom and top quarks, such as the top assisted technicolor 
model, but it is beyond the scope of our review. 

Since charm quark is indeed sufficiently heavy (will discussed below, "heavy" means it 
is heavier than the binding energy scale Aqqd), it is natural to use the potential model to 
evaluate the spectra of J/ip and other family members: r]c, Xc and even he etc. and their 
excited states, as well as their fine-structures and decay modes. However, since it is not 
too heavy, the relativistic correction to the potential which generally consists of a Coulomb 
piece and a confinement one (for example the linear potential) , is more serious than that 
for T family. Even though, by adjusting parameters, great success has been achieved for the 
heavy charm-quarkonia (charmonia) and the results are satisfactorily consistent with data 
for ground and lower excited states. Very recently Voloshin gives an enlightening review 
on charmonium where relevant topics are discussed in details 7|. However, the story is far 
from its end yet. There have been many puzzles in the field, especially as more accurate 
measurements are done, they do not disappear. Moreover, several new resonances have been 
observed and they seem not easy to be described by the simplest valence quark structure, 
i.e. meson is composed of quark-antiquark and baryon is composed of three quarks. 

Now, let us present a rough list about the puzzles. The first one may be the famous 
pvr puzzle where the branching ratio of tp' pir is too small compared with 8|. Then the 
decay mode J/tp ^ pn is forbidden by the hadronic helicity conservation, thus its rate must 
be sufficiently small, but by contraries, it is one of the main decay channels of J/ip. The 
sizable D'^ — mixing which was recently observed j^, indicates that there must exist a 
mechanism beyond the Standard Model (SM). There are many newly observed resonances, 
which may demand interpretation. 

On the theoretical aspect, great efforts have been done to look for reasonable explanations. 
The first step is to make sure whether the puzzles in the new observations are due to the flaws 
of our theoretical framework or there is new physics beyond the SM. Indeed, to understand 
the experimental measurements, one must calculate the corresponding quantities, where 
evaluation of the hadronic matrix elements is the key and obstacle. Since the hadronic 
matrix elements are fully governed by the non-perturbative QCD effects which at present 
there is not a reliable way to deal with, all the results we have achieved must possess certain 
uncertainties. To properly estimate the errors and reliability of the results becomes an 
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important issue in theoretical calculations. 

Below, we will discuss all the topics in separate sections. We first recall the mechanisms 
which govern the weak transitions of D mesons, and then focus on discussion of the rare 
decays because the regular Cabibbo-favored channels have been thoroughly investigated in 
both theory and experiment. Then since all the topics are related to theoretical evaluation 
of the hadronic matrix elements, in next section we review the status. In section llVt we 
concern the final state interaction, especially focusing on the hadronic triangle calculation 
and briefly discuss other schemes. In section |Vl we will discuss the hadronic helicity selection 
rule and its violation. In section |Vl] we concern the pvr puzzle and in section IVIIl we show 
how the QCD multi-expansion theory works very well for the pion radiation from excited 
states of T, but has difficulties for il){nS) il){mS) + tttt (n > m). In section IVlIU in 
some details, we discuss the newly observed resonances at the charm energy region and 
some of them seem to be exotic and need a reasonable interpretation. In section IIXI we 
discuss the — mixing and related theoretical proposals, meanwhile we also concern 
the possible CP violation observable. In section |Xl we consider the charmed baryons and 
double heavy baryons Hcc, where the difference between the lifetimes of Ac and D"^^ is 
especially concerned and we will also briefly discuss the charmed pentaquark. In section 
Kit we specially discuss the diquark structure in baryons, especially in the heavy baryons, 
because it is an important issue for studying baryons. The last section is devoted to a brief 
discussion. 

II. THE EFFECTIVE LAGRANGIAN AND TRANSITION AMPLITUDES 

Generally, this is a topic which is familiar to most of theorists working in this field. Thus 
after a short introduction, we will turn our attention to the application of the theoretical 
framework for investigating rare decays of D mesons and J/ip. 

A. The decay constants of Ds and D^'^ 

The most important parameters for the weak decays of D and Ds mesons are their decay 
constants. Generally the decay constant of a pseudoscalar meson is measured via its leptonic 
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decays P ^ lu (l = e, fx, t as long as kinematics allows.) and the width is written as [K 



r(P ^ Iv) = ^flm]Mp{l - ^)\V,,.\\ 



where Mp,mi are the masses of the pseudoscalar meson and the lepton, Vqq' is the corre- 
sponding CKM entry and q, q' are the valence quarks in the pseudoscalar meson. There fp 
is the decay constant which we are going to obtain. The CLEO collaboration has achieved 

fD+ = (222.6 ± I6.7tli) MeV. 



in the decay of ^ 11 1 
Rosner 



10|] has obtained the average of fp,^ as 

fp^ = (274 ± 10) MeV. 
There is a discrepancy of about three standard deviations between this result and a recent 

n 

unquenched lattice calculation [14l |. 

Very recently, the CLEO collaboration reported their new result as [li 

fD+ = (205.8 ± 8.5 ± 2.5) MeV, 

by assuming \Vcd\ = \Vus\- They also obtained 

B{D+ n+v) - B{D- n-v) 



0.08 ±0.08. 



B{D+ ji+u) + B{D" ii-p) 
which means that no CP violation was observed in the leptonic decay. The Belle group 



reported 



12| 



/^^ = (275 ± 16(stat) ± 12(syst)) MeV. 



Only when the decay constants are accurately measured, the theoretical predictions on 
the hadronic transitions can be trustworthy, so that more precise experiments are necessary. 
In fact, measurement on decay constants is not easy, not only because the so far available 
database on especially Ds is not large enough to guarantee a high statistics, but also 
reactions, such as P ^ /i^ + 7 can influence extraction of fp. The BESIII will collect the 
largest database of D mesons, so one may expect to get very accurate decay constants of D 
and Do. 
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B. The effective Lagrangian of weak interaction 



The charmonia J rjc and their excited states mainly decay via strong interaction which 
is the OZI suppressed (will be discussed in later sections). Instead, D mesons decay via weak 
interactions or electromagnetic radiation. In this framework, the semileptonic decays are 
well understood, so that we concentrate ourselves here on the non-leptonic decays. Let us 
briefly review the general situation of the D decays. The effective Lagrangian for weak 
interaction is written as 



-|Ac|=l 



G 
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where q can be either d or s, and the flrst operator Oi 



i=3 



+ h.c, 



{sc)v-A{uq)v-~A with {qq')v-A = 
g7^(l — 75)^' originates from the tree level while others are induced by loops, O3 to Oq are 
the strong penguin operators and the rest = 7 to 10) are due to the 7, Z— penguins and 
box diagram 15|. 



C. Rare decays of D meson 

The weak decays of D mesons can be categorized into Cabibbo favored, Cabibbo 
suppressed and doubly suppressed modes. The flrst type includes the modes such as 



K'^ + TT'^ etc. and the second one was discussed by Abbott, Sikivie and Wise 16|. For 
the Cabibbo favored decay modes, one usually only considers the so called spectator mech- 
anism where the light quark behaves as a spectator when the heavy charm quark transits 
into s quark plus a ud pair. In this picture annihilation and VT-exchange between charm 
and light anti-quark can be neglected because of the linear momentum matching (namely, 
for annihilation and VT-exchange, a quark-anti-quark pair must be produced from vacuum, 
or in other words are produced by soft gluons, so that their linear moneta are small whereas 
the quark- antiquark pair occurring directly from the effective vertex possess large linear 
momenta. As they combine a quark (anti-quark) emerged from vacuum to constitute a 
hadron where the two constituents must have close linear momenta the large momentum 
difference would greatly suppress the probability.). However, for the Cabibbo-suppressed or 
even doubly-suppressed modes, the annihilation and IV-exchange mechanisms and as well as 
the penguin contributions become important. In fact, there are several channels where the 
spectator mechanisms do not contribute at all, thus these modes would be ideal places to 
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study such small effects which may manifest some unknown mechanisms. Moreover, for the 
CP violation, the contribution from the penguin and even the electro-weak penguin would 
be crucially significant. We will discuss these issues in the following sections. 

Another interesting rare decay mode are those processes, where the light quark (anti- 
quark) transits while the heavy charm behaves as a spectator which only provides a color 
source. Such reaction includes D* ^ D + 7, D* Dn etc. where the photon and pion 
can be emitted from either charm or light flavor. We used to study a special case that a 
heavy baryon containing two heavy quarks radiates a photon and transits into a lower states 



with the same flavor 



171 ] in terms of the Bether-Salpeter (B-S) equation. It was indicated 



that the branching ratios of such decays are very small and hard to measure at the present 
luminosity, however, for BESIII and LHCb, the situation may be greatly improved. 

Recently, Li and Yang calculated the branching ratios of + + u, — »• 

+ + u, D'^ + + e~ in SM, however, their results indicate that only the 

branching ratio of Df + + u could reach lO"*^. According to the sensitivity of 

BESIII, B-factories, Super-B and LHCb, it might be observed at Super-B and LHCb, but 
not at others. On other side the observation may offer a probe for testing the working 
mechanisms which govern the behaviors of the light flavors in hadrons which are usually 
treated as passive spectators in most of reactions, as aforementioned. 

An interestin g d iscovery draws attention of theorists, it is the observation of baryonic 
decay Dg pn ly], which can only occur through the ly-annihilation topology, so it was 
supposed to be very suppressed as aforementioned for the meson case [2^. Chen et al. 21] 
indicated that the short- distance contribution can only make the branching ratio as large 
as 10~^ which is much smaller than the data. Thus they suggested that the long-distance 
contribution via the FSI can enhance this value, so that they claimed that it is a dynamical 
enhancement of the VT-annihilation topology in Dg decays. It is worth further studies indeed. 



D. Weak decays of J/ip 

This is another type of rare decays which may provide us with some information about 
the structure of J/tp. Generally J/t/,' would decay via strong or electromagnetic interactions. 
The strong decay is realized via a process where the constituents c and c annihilate into three 
gluons which eventually fragment into hadrons, it is an OZI-suppressed reaction and that is 
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also why J/ip is a narrow resonance and evaded observation before 1974. It is believed that 
the decay width is proportional to the wavefunction of J/tp at origin which can be easily 
obtained by measuring its leptonic decay width. However (see below), a violation of the 
hadronic helicity selection rule indicates that such a picture may be not completely correct, 
therefore to investigate the structure of J/ip (if it has a hybrid component etc.), study on 
the weak decay of J/ ip might be very helpful. We calculated the branching ratio of the semi- 



leptonic decay J ftp +e+ + |22| in the QCD sum rules and obtained it to be of order 

of 10"^'^. Then with the gained parameters, we extended our calculation to the non-leptonic 
weak decays of J/ip. The results show that the branching ratio of inclusive weak decays can 
reach order of 10~^, and a special channel J /ip D^*^ + p 23| has a larger branching ratio 
of about 5.3 x 10~^ which might be measured by BESIII, B-factories, Super-B and LHCb, 
as we wish. 

In fact, all such rare decays may be important for better understanding of the structure 
of J/ip and the governing dynamical mechanisms, even though accurate measurements on 
them are extremely difficult. We lay our hope on the very large database of the facilities 
which will be available soon. 



III. HADRONIC MATRIX ELEMENTS 

This is probably the most difficult problem in hadron physics which almost covers the 
whole field of high energy physics, and is definitely confronted by anybody. The reason 
is that hadronization occurs at the energy scale below Aqcd where non-perturbative QCD 
effects dominate and so far there is no an effective way to accurately evaluate the effects 
yet. Much efforts have been made to handle the problem. The simplest way is using the 
naive factorization where a hadron, generally a meson, is emitted and can be factorized out 
from the hadronic transition of one hadron (meson or baryon) to another one. Then the 
rest transition amplitude can be parametrized by a few form factors which are obtained by 



fitting data [2J]. The transition matrix element can be analytically decomposed into a few 
terms according to the Lorentz structure and the parity conservation because hadronization 
is a process where only strong interaction applies. The advantage of this method is that 
it is simple and since the form factors are fixed by fitting data, they can be extensively 
applied to study the weak decay rates. It is simple and consistent with data within a rather 
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wide range, however, there are obvious shortcomings. First the factorization is not always 
legitimate, as Buras et al. pointed out 25|, the matrix element of operator 



cannot simply written as the factorized form because a term proportional to 

would appear and phenomenologically it causes an effective N^. and the coefficients are 
deformed as 

Ci + C^/Nlff or Ci + C^/Nl^f. 



where N^^f is no longer 3 26 1. 



Moreover, such factorization is based on the spectator mechanism where the transition 
occurs at the heavier quark leg and another light component would play a role of a spectator. 
In this way, the annihilation and H^-exchange sub-processes are not properly included. On 
other aspect, as we know, the annihilation and H^-exchange sub-processes might be impor- 
tant, when dealing with the inclusive processes, especially for evaluating the lifetimes of D 



mesons 271 ]. 



Thanks to the heavy quark effective theory (HQET) where an extra symmetry SUf(2) x 
SUs{2) is considered, one can reasonably evaluate the transition between two heavy mesons 
containing 6 or c quarks. It has already becomes a criterion for testifying validity of any 
theoretical calculations where hadronic matrix elements are evaluated, as if the heavy quark 
limit is taken, the results must be qualitatively consistent with that obtained by the HQET. 
Moreover, Georgi generalized this scenario for dealing with transition between two h eavy 



baryons, each of which contains two heavy quarks in terms of the superflavor symmetry [28] . 
On another side, charm is not heavy enough to be treated as a real heavy quark, and the 
l/rric corrections may be important, even for heavier b-quark, the 1/mf, corrections are not 



negligible in practical computations 29|]. Therefore, just as the results under laige-Nc limit 
in the 1/Nc expansion theory only possess qualitative meaning, the obtained values under 
heavy quark limit (i.e. let mq oo) correspond to the leading order, corrections must be 
accounted while comparing with more accurate experimental data. There are many works 



to consider how to properly evaluate the l/mg corrections 29] 
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In the recent years, some effective theories have been developed to justify the factoriza- 
tion, especially in B physics. The perturbative QCD (pQCD) is based on the factorization 
theorem. It states that for a process with large momentum transfer, the physical amplitude 
can be factorized as a product (f)®H where H is the factor corresponding to the quark-level 
hard process amplitude which can be calculated in perturbation theory order by order and 
stands for the soft part. The later one is not calculable in the perturbative way. In general 
the product is related to a convolution integration over the wave functions of the initial 
and final hadrons. In the literatures, there are two different versions of pQCD to deal with 
the factorization, one is the familiar coUinear factorization and another is the so-called kx 
factorization. The crucial difference between the two schemes may be that factorization 
is possible to treat the problem of endpoint divergence. The details of the two schemes can 
be found in relevant literature 30|. The soft-coUinear effective theory (SCET) is a recently 



developed effective field theory to simplify the processes containing light energetic hadrons 



31 



32[. The great development is that the proof of the factorization theorem can be per- 
formed at the operator level in the SCET which is much simpler than the diagrammatic 
analysis in pQCD. However, on the other hand, except very few processes, such as the pion 
transition form factor, have been proved to be factorizable, most exclusive processes can- 
not be rigorously proved to be factorizable. Moreover, the factorization proofs are usually 
limited to the leading order and therefore the application factorization theorem is still an 
assumption. The factorization may be applicable for the decays of bottomed mesons and 
baryons, but for charmed mesons D and baryons, it is indeed questionable. In D meson 
decays, the energy of final light meson is at the order of Aqcd which is not high enough to 
perform a perturbative analysis. Another question arises from the substantial corrections 
in power of A^qcd /rric. All these facts make it difficult to apply pQCD or SCET into charm 
decays. In another approach named as the Transverse Momentum Distribution (TMD), 
where factorization is performed by taking transverse momenta of partons into account. For 
example in the reaction tt*^ + 7* — > 7 the form factor can be written as a convolution integral 

F(g') ~0®^®i^(i + c(g-'). 

where is the light-cone wavefunction of the pion, S is an additional soft factor and H 



corresponds to the hard scattering part 33[|. 
A very recent work by several authors 



34( 1 indicates that for the kx factorization scheme, 
10 



at loop-level an extra term which is related to the so-called light-cone divergence appears 
when a unitary gauge is employed in the calculation and disappears in the Feynman gauge. 
This explicitly manifests that the the /cy factorization at loop-level is not gauge-independent, 
so that is violated. This statement is still in dispute. Even though the kx factorization 
cannot be a strict theory according to the field theory, it can definitely treated as a successful 
phenomenological model and is applied to calculate the transition amplitudes where heavy 
hadrons are involved. Therefore generally we can trust the theoretical results which are 
obtained based on the /c^ factorization. Very, very recently, H. Li argues that the work 34 1 
might make some calculation mistake and he presented a result which is free of the light-cone 
singularity, so the gauge invariance of the kx factorization is kept [35|. Since it is a serious 
dispute, we will follow the further development in the interesting regime. 

Besides these methods which may stem from the quantum field theory, there are some 
traditional methods which have been widely applied to calculate the hadronic matrix el- 
ements including for example: the harmonic oscillator model 36||, the constituent quark 
model [s^], the constitute quark meson model (CQM) 38, SoJ, light front quark model 40|. 



color-singlet model 41[, color-octet, and color evaporation 42], especially the non-relativistic 
QCD (NRQCD) where an expansion in powers of velocity of the heavy quark v is naturally 
made 43|. Besides these phenomenological models whose parameters must be fixed by fit- 
ting data, the theoretical framework QCD sum rules [44 1 is based on quantum field theory 
where only the perturbative vacuum is replaced by the physical vacuum. Because of the 
properties of the vacuum, a series of condensates of quark-pair, gluons and quark-gluon etc. 
are introduced to describe the non-perturbative QCD effects. It has achieved great success 
in phenomenology. However, on the other side, it is an extrapolation from the region where 



perturbative QCD works reliably 



451 ]. In the expansion only the operators with lower dimen- 



sions are retained and moreover, to extract physical results a reasonable plateau is required, 
where the threshold values are determined, thus an error of about 15% is unavoidable. 

Quite amount of phenomenological models have also been employed to calculate the 
hadronic matrix elements whose energy scale is Aqcd- As well known, the hadronization 
is fully determined by the non-perturbative QCD effects and so far there is no a reliable 
way to evaluate them based on quantum field theory or any other first principles. One 
of the goals of our research in fact, are to determine the mechanism, which governs the 
reaction, generally the fundamental theory is the standard model (SM) which at present 
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no one doubts due to its remarkable success, thus we can reliably (if ignore contributions 
from new physics beyond the SM) determine the hard factor due to the asymptotic freedom 
of QCD. To extract important information, such as determining the CKM matrix elements 
and checking the unitary triangle, exploring CP violation, one indeed needs to have a more 
accurate estimation of the hadronic matrix elements, otherwise the physical picture would 
be contaminated by the inaccuracy. Therefore, to understand the physical world, reliable 
estimate on the non-perturbative effects are absolutely necessary and all the efforts along 
the line are worthwhile. 



Moreover, one can also use the MIT bag model with ta^ 



dng into account the recoil effects 



47( 1 to describe the wavefunctions 



46(1, the chiral bag model and even the fiux-tube model 
of the initial and final hadron states when carry out the calculation of the transition matrix 
elements. 

Indeed all the models have their own reasonability and advantage, but there are obvious 
flaws and un- reasonability, and because they are not coming from a basic principle, one 
can never expect that they can be perfect. Therefore on one side, even though the models 
are not perfect, they have applicability and if they are properly applied, reasonable results 
should be reached. 

As a conclusion, estimation of the hadronic matrix elements is crucially important, but 
so far, there does not exist a way to fulfil the job yet and one can only apply the available 
models to estimate them with certain reliability. 

IV. FINAL STATE INTERACTION 

Besides the estimate on the hadronic matrix elements which are directly related to the 
transition, there are secondary reactions, namely the final state interactions (FSIs) which are 
also very significant for charm-hadron decays. Such processes are due to strong interaction 
and occur at hadron level, thus also cannot be derived by perturbative QCD. Fortunately, 
one can use the chiral Lagrangian to evaluate the long-distance effects and we will discuss 
this issue in this section. Some phenomenological models had been suggested to estimate 
;he FSIs in D meson decays: one-particle-exchange model 48|] and the Regge pole model 



49|. 



The final state interaction (FSI) in the charm-tau energy region is very important 50, iM] 
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According to the concept, the FSI can be categorized into the quark-level and hadron level 
FSI processes. The quark level FSI process refers to the quark interference due to the 
identical fermion statistics. It was noticed by Stech et al. long time ago to explain the 
lifetime difference of and and we will come to this subject in later sections. Now let 
us concentrate on the second category of FSI i.e. that at hadron level. In those processes, 
the initial hadron first decay into intermediate hadrons (usually two hadrons) and the two 
hadrons would re-scatter into the final states. Since the re-scattering occurs via strong 
interaction, the isospin must be conserved. 

The re-scattering occurs at hadron level and both of the hadrons are in color- singlet, 
thus the interaction between the hadrons cannot be described by one or even a few gluon- 
exchange and it makes the whole calculation more tricky and uncertain. The responsible 
effective theory in this field should be the chiral Lagrangian. However all the coefficients in 
the lagrangian cannot be obtained from an underlying theory, such as QCD at present, and 
can only be fixed by fitting data, this brings up very serious problem and uncertainties for 
the theoretical evaluation. 

To evaluate the re-scattering effects, some authors suggested to calculate the absorptive 
part of the triangle diagrams whose internal legs are intermediate hadrons and external 
lines correspond to the initial hadron and the final daughter hadrons 



In fact the 



absorptive part of the triangle corresponds to the real Final State re-scattering because 
the two intermediate hadrons which are directly coming from the initial decaying hadron 
are on their mass shells. The hadron(s) exchanged between the two intermediate hadrons 
at t-channel not only needs to possess proper quantum numbers, but also has to conserve 
energy-momentum, so that it is obviously off-shell. As at the triangle apexes we apply the 
effective interaction vertices which are extracted from the chiral Lagrangian, one needs to 
introduce phenomenological form factors to compensate the off-shell effects. Usually there 
are various types of the form factors which are widely adopted in literature, the simplest 
one is the pole form^. 



A^ — m 

A2-g2 



2 



The usually adopted form of the form factor is )" where n = 1 is the monopole form, n = 2 is the 

dipole form and as n > 2, it is a multi-pole form which is seldom selected in literature. Besides, there are 
exponential and other forms. 



13 



where q and m are the momentum and mass of the t-channel exchanged hadron and A is a 
phenomenological parameter whose value is beheved to be close to 1 GeV. It is noted that the 
form factor can also paly the role of the cut-off in the Pauli- Villas renormalization scheme, 
so that in the calculations, no ultraviolet and infrared divergences bother us. On other 
aspect, this also leads to parameter-dependence and makes theoretical predictions uncertain. 
Therefore, this calculation can only tell us the order of magnitude for the concerned reaction 
unless we can use some data as inputs to fix the model parameters. 

The importance of FSI can be understood as one studies the decays of K^K^ 



and K^K~ j49|], the former process can only occur via ly-exchange diagram which 

is very suppressed according to the general analysis, moreover, due to a CKM cancelation, 
this reaction should be proportional to an SU(3) violation, so that should be very small, 
in comparison, the later one is a Cabibbo favored external emission process and should 
be overwhelmingly larger than the former one. However the data show that B{D^ — >■ 
K+K-) = (3.84 ± 0.3) X lO^^ and B{D° 2A'°) ~ (3.7 ± 0.7) x 10-^ which implies 
B{D^ —>■ K^K^) is comparable with B{D^ —>■ K^K~). This can be easily realized via a 



re-scattering K^K — s> K K . In our work (49|], we showed that it can be realized with the 



data measured in experiment on KK scattering as inputs. This simple example confirms 
the importance of FSI in charm physics. 

Moreover, the final state interaction provides a strong phase which may lead to CP 
violation. As well known the direct CP violation in decays is 

T{A ^ B)- r{A 5) ~ sin(ai - aa) sin(0i - ^a)- 

where B may be a CP eigenstate and a, (j) are strong phase and weak phase respectively. 
It is obvious that there at least exist two independent channels which have different strong 
and weak phases, otherwise the direct CP asymmetry is zero. Due to the final products 
may originate from another weak process and therefore the reaction can have different weak 
phase from the direct transition and the FSI can provide a strong phase, i.e. the phase shift 
in the language of scattering, thus their interference can result in a non-zero CP asymmetry 



|5l|. 



To calculate effects of the final state interaction, there are two possible ways to adopt, 
i.e. the Regge-pole model and the hadronic loop. It was discussed that the Regge model 
may apply in higher energy region whereas the hadronic loop method is more suitable in 
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lower energy region. To certain accuracy they are consistent. 

Generally, as many authors discussed, the two intermediate hadrons are real on-shell 
particles, therefore, so that one only needs to calculate the absorptive (i.e imaginary) part 
of the triangle diagrams. In fact, it is the real final state interaction by the common sense. 
However, as Suzuki jsS] pointed out, the dispersive part of the triangle can also paly a role 
to influence the transition amplitude. By calculating the dispersive part of the triangle, 
we evaluate the branching ratio of J/tp ^ PV where P and V stand for pseudoscalar and 



vector mesons 



54j ]. This is related to the famous pvr puzzle which will be discussed in later 



sections. Our strategy is that we calculate the dispersive part of the triangle by keeping the 
parameter A as a free parameter to be determined, and by assuming the SU(3) symmetry we 
set the direct transition rate as another free parameter, then by fitting two special channels 
{J/ip P^T^^i K^*K~ which are more accurately measured) we obtain the two parameters. 
With them we calculate the branching ratios of other channels and find that the results 
which are listed in the table below, are well consistent with data. 



Dccav mode 



BRxlO ''(Expcrimcnt1f571 



eg'' (10- ^ GcV- 1 ) (Theory) 



g-f'^(10-^GcV-l)(Thcory) 



BRxlO ^(Theory) 



ip-q 



4>V 



0.65 ± 0.07 0.33 ± 0.04 1.58 ± 0.16 0.167 ±0.025 



2.08 ± 0.25 1.65 ± 0.26 0.89 ± 0.0 



1.27 ± 0.13 0.46 ±0.069 



6.44 



6.01 



3.47 



8(fitting) 1.65(fitting) 0.93 0.61 



4.2(fltting) 5.0(fltting) 0.71 



TABLE I: The first two modes are well measured and the theoretical model parameters are obtained 
by fitting them [5j]. 



It is interesting to note that usually only at lower energy region, the FSI effects are more 
significant, but sometimes, the small effects may be also non- negligible. This is the case for 
studying CP violation. As indicated, the direct CP violation is induced by an interference 
between at least two channels which have different weak and strong phases. Even though the 
FSI in certain cases is much smaller than the main contribution which usually comes from 
the tree level and does not possess a strong phase, FSI then provides a non-zero strong phase 
and its weak phase may be completely different from the tree level one, thus an interference 
between it and the tree amplitude would result in a non- zero CP violation. We calculated 



such a possibility for the Be decay 



and conclude that in the case the pQCD calculation 



is still valid and the FSI effect only raises a minor contribution to the total decay amplitude. 
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but its strong phase is not zero, so may help to build up an observable CP violation effect 
as it interferes with the tree contribution which is calculated in pQCD. 

This indicates that the FSI is important, but as aforementioned, the large uncertainties 
of the whole scenario and parameters involved in the calculations make the theoretical 
prediction not accurate, and all these are worth further studies. To make the theory in a 
better shape, we need more information from experiment. 



V. HADRONIC HELICITY SUPPRESSION IN J/0 DECAYS 

It is generally believed that the narrowness of J/ip is due to the Okubo-Zweig-Iizuka 
(OZI) rule 561]. Violation of the OZI rule would give rise to a non-zero transition rate. 
The main two-body decays of J/il) can be categorized as J/il) PP', J/tp ^ PV and 
J/tp —>■ VV where P, V stand for pseudoscalar and vector mesons. Besides, definitely there 
are other modes, for example, the scalar and axial vector final states, and as well, there 
is also possibility to decay into two baryons such as pp, AA etc. but the corresponding 
branching ratios of such channels are rather small and the result is well understood in 
our theory. To get a better understanding of the OZI rule, people turn to study the radiative 



decays of orthoquarkonia [58|] where only a hadronic transition matrix element is needed. 
With certain approximations they obtained numerical results which were roughly consistent 
with data. The process involves a five-point Green's function and the Feynman integration 
is complicated. Thanks to the developments of the calculating techniques, we carried out 



a full calculation 
authors of Ref. 



59| and the results are qualitatively consistent with that obtained by the 
58l | and the data available at that moment. The success indicates that our 
knowledge on the OZI rule might be right with a tolerable error. It is also noticed that in 
the radiative decays, the hadronic effects of J ftp are included in its wavefunction at origin 
which is obtained by fitting data of the leptonic decays of J/ip and it must be accurate 
enough. 

To further testify the OZI rule, we calculate the process oi J/ip ^ tttt which is an isospin 
violating one. Usually it is supposed that this process is induced by the electromagnetic 
interaction, i.e. via J/t/.' ^ 7* ^ tttt. However, as well known, there are two sources for 
isospin violation, one is the electromagnetic interaction and another is the mass splitting 
of u and d quarks. Considering the mass difference of u and d quarks, we calculate the 
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OZI process J/ip ggg mr 60j and got its amplitude which is comparable with the 



contribution from electromagnetic interaction. Concretely we find the OZI amplitude is 
proportional to (m„ — md)/Mj/^ which clearly manifest the isospin violation. We should 
will carry out similar calculation oi J/tp ^ '-/* ^ vrvr in the same scenario and compare it 
with the OZI process in our later work. Then we continue to calculate the branching ratio 
oi J/ip pn which is supposed to be fully dominated by the OZI process J/ip ^ ggg pn. 
We obtained a branching ratio which is one order smaller than the data More discussion will 
be presented in the next section). 



In fact, it is not a surprise, because a long time ago, Brodsky and Lepage [61[ indicated 
that as the vector-gluon coupling conserves the quark helicities, the total hadronic helicity 
is conserved and can only be violated at order of m/Q or higher, where m is the light 
quark mass in the final hadrons and Q is the transferred momentum scale. Our numerical 
results are consistent with this rule, i.e. the process is suppressed by the hadronic helicity 
conservation and we explicitly show that the amplitude of the transition is proportional to 
{niu + nidj/Mj/^ which confirms the observation of Lepage et al. However, the allegation 
sharply contradicts to the data where the branching ratio oi J/ip pir is (L69 ± 0.15)% 



571] and is one of the dominant hadronic modes. To compromise this obvious discrepancy. 



Branbilla 



6l| suggested that either J/ ip contains other constituents, for example is a hybrid. 



or there is some unknown mechanism which violates the hadronic helicity conservation and 
results in a larger transition amplitude. However both of the interpretations would receive 
very crucial challenges because for a long time, people believe that J/ip is composed of 
charm and anti-charm quarks and almost overwhelming works are based on this picture. 
The second one also does not seem to work, because if it were true, the mechanism would 
exist in other channels and infiuence all theoretical predictions. 

Associating with the phenomenon where ip' — ^ pn is peculiarly suppressed, we are inclined 
to the first interpretation, namely J/ip is not a pure c — c bound state, but ip' is. We will give 
more discussions in the later section about the pir puzzle. On other aspect we used to follow 
Suzuki and consider the long-distance contribution to J/ip —* pir and our results show that 
it is possible to explain data, but the answer is still not satisfactory yet (see later section 
for more discussions). We used to follow Suzuki and consider the long-distance contribution 
to J/tp ^ pn and our results show that it is possible to explain data, but the answer is still 
not satisfactory yet (see later section for more discussions). 
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VI. THE pvr PUZZLE 



As widely discussed, the puzzle has been raised for a long while. In the regular theoretical 
framework, there should be a relation 

BRji^'^ggg) _ r(^- ^ e+e") r,(J/V;) 
BR{J/ij^ggg) r(J/^ ^ e+e") Tt{^') ' 

where Ff is the total width. This ratio comes from the fact that if both J/tp and tp' are 
c — c bound states, in the hadronic decays, c and c annihilate into three gluons which then 
convert into hadrons, whereas in the leptonic decays, c and c annihilate into a virtual photon 
which turns into a lepton-pair. In this picture, the amplitudes of the hadronic decay which 
occurs via a three-gluon intermediate state, and the leptonic decay which occurs via a virtual 
photon intermediate state are proportional to the wavefunction at origin iplO). If everything 
worked well, the ratio should be close to 12~ 14%, which is called as the 14% rule (now, 
it is sometimes called as 12% rule, anyhow it is a sizable number in contrast to the data.). 
However the data tell us that this ratio is much smaller than this value. 



Some theoretical interpretations have been proposed. Rosner et al. [62] suggested that 
the quantum number of the observed ip' may be not a pure 2S state which is the first radial 
excited state of the cc system, but a mixture of 2S and ID states. The amplitudes are 
instead 

{pn\ij') = {pn\2^Si) COS0 - {pn\l^Di) sin0 ~ 0, 

(pvTl^/'") = (p7r|235i)sin0+ (p7r|l^Di)cos0~ (p7r|2^Si)/ sin 0, 

where is fixed as —27° or 12° by fitting data. By the destructive interference between 
the contributions of the two components to the amplitude of ip' — pvr, the smallness is 



explained. Suzuki [53|] alternatively suggested that the relative phase between the one- 
photon and gluonic decay amplitudes or hadronic excess in ip' decay may result in the 
small branching ratio. The final state interactions may also give a reasonable explanation 



63( 1 . The first proposal can be tested in the decays of tp" pvr which has not been 



well measured yet. In Ref. (53|, the author suggested that the one-photon amplitude is 
sizable and it can be tested in some other modes, for example vrvr if the process is 

dominated by the electromagnetic interaction. The hadronic excess can also receive tests in 
the decays of other higher excited states of ip— family and even T— family. The final state 
interaction may play an important role in D and B decays, and also in decays of tjj mesons 
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as suggested in literature. The difficulties are how to properly evaluate such effects. The 
final state interaction process is induced by strong interaction at lower energy region, thus 
it is governed by the non-perturbative QCD which is not fully understood in the present 
theoretical framework yet. People need to invoke some phenomeno logical models to carry 
out the calculations. We will give a more detailed discussion on estimation of final state 
interaction, here we only use our results to discuss the puzzle. 



In our work 



5J], we simultaneously consider the FSI and the direct decay of J/ip into a 



vector and a pseudoscalar mesons and conclude that both of them contribute to the widths 
and their interference should be destructive to explain data. This observation indicates that 
even though the OZI-forbidden process is sizable, it cannot be consistent with data. The 
result implies that the hadronic helicity conservation indeed greatly suppresses the process 
of J ftp ^ pn and as the data demand an explanation, one should consider what is the origin 
of the problem. 

In a straightforward calculation based on the SM, we estimate the decay width of the 
OZI forbidden process J/tp ^ pn [60[], and find that the width is indeed proportional to 
{Triq/mji^Y which is coming from the hadronic helicity suppression factor. Numerically the 
branching ratio of J/tp ^ pn should be smaller than 0.1%. The same situation appears 
for ip' pn. It was qualitatively discussed by Brodsky et al. 6l|. As aforementioned 
in last section, to testify the calculation, we recalculate the subprocess J/tp ^ 3g ^ nn, 
which is an isospin violating reaction and usually is supposed to be dominated by the 
subprocess J/ip ^nn because the EM interaction violates isospin as well known. Our 
result indicates that in the OZI forbidden subprocess the transition amplitude is proportional 
to (m„ — Tn(i)/mji^, i.e. the mass difference results in the isospin violation instead. Our 
numerical result is of the same order as the contribution from J /ip '-/nn. All the results 
are consistent with our physics picture and qualitatively reasonable. Therefore we can trust 
our calculations for the process J/ip —* pn. Our numerical results are listed in Table HTl 

As indicated in jfill , the structure of J/ ip may be not a pure cc charmonium, but consists 
of other components, such as hybrids ccg, ccqq and etc. 

To understand the smallness of the ratio R, one can expect that either there is a problem 
with ip' as Rosner et al. do, or something obscure in J/ip structure as Brodsky and many 
others indicated. Our above numerical results show that even though the hadronic selection 
rule works in the cases of J/tp ^ pn and tp' ^ pn, the suppression is not too serious and 
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TABLE II: Decay widths (T) of J/tlj ir'^p + n based on the three distribution functions, 



b2 and 03, respectively 



muiMeY) 


md(MeV) 


r(<^i)(MeV) 


r(02)(MeV) 


r((/>3)(MeV) 


exp(MeV) 


2 


2 


1.04 X 10^"^ 


7.21 X 10^^ 


5.11 X 10-"^ 




3 


3 


2.36 X 10-"^ 


1.6 X 10-^ 


1.17 X 10^3 




4 


4 


4.12 X 10"'' 


2.9 X 10"^ 


2.08 X 10-3 


(1.06 ±0.08) X 10-3 


5 


5 


6.69 X 10^"^ 


4.54 X 10^"^ 


3.38 X 10^3 




6 


6 


9.75 X 10-"^ 


6.68 X 10""^ 


4.88 X 10-3 





the theoretical prediction is only one order smaller than the data. 

Therefore a tentative conclusion may be drawn that the pvr puzzle may be not due to the 
mixing structures of ip' and ip" , but neither to an anormal structure of J/ ip itself. It seems 
that both of the proposals cannot independently explain the "puzzle", more complicated 
mechanisms may be needed. 

This is a great challenge to our understanding because the cc structure of J/ip has been 
recognized almost from very beginning of its discovery. If it is not a pure cc, all the previous 
works in terms of the potential models where many parameters are fixed by fitting data 
should be re-considered. Or there may be some other mechanisms which were not taken 
into account, or may exist contributions from new physics beyond the standard model (SM). 
However, the later seems not very promising because the concerned energy range is rather 
low and SM works perfectly well to explain the data for most states and processes. Thus it 
is obviously inclined to the first proposal that J/ip is not a pure S-wave bound state of cc. 

It is noted that the calculated results unless for the distribution 03, are one order smaller 
than the data. The same situation happens to the but is still hard to draw a conclusion 
that the 14% rule is due to existence of higher Fock states in J/ or other mechanisms which 
further suppress the reaction of ip' —* pir, may be both. It forms an intriguing challenge to 
our understanding of the hadron structures. This whole picture also applies to T family, 
therefore the future experiments would provide hints to finally solve the puzzle. 

There have been some theoretical explanations besides that we discussed above, in the 



work by Mo, Yuan and Wang 



641 ]. the authors described the recent status of theoretical 
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research as well as the experimental measurements on the interesting subject. 

Interesting, there is also an alternative opinion towards the subject, Suzuki 
65l | deny it as a "puzzle", because they considers that the electromagnetic interaction may 



531 ] ■ Zhao 



play an important role in J/ tp decays where cc annihilate into a virtual photon which later 
fragment into hadrons. In the picture, it is supposed that a destructive interference between 
the contribution of three-gluon and single-photon processes would suppress ip' pvr. Since 
its amplitude can be roughly estimated in terms of the measured rate of J/ ip leptonic decay, 
there should be a strong constraint on the proposal. Moreover, if it is true, the interference 
would also appear to other decay modes and the picture should be further investigated 
and tested by more accurate experimental data available in the future, especially from the 
BESIII. 



VII. QCD MULTI-EXPANSION AND STRONG RADIATIONS 



In an enlightening paper by Kuang 66|, in details, the author explained the work of 
Yan and Kuang [67| where they successfully initiated and developed a complete theoretical 
framework, the multi-expansion method in QCD. The theory properly deals with emission 
of light hadrons during heavy quarkonia transitions. Concretely, one mainly studies the 
processes such as T{nS) T[mS) + vr + vr or ipijiS) — > tp{mS) + tt + vr with n > m. In 



Ref. 66[, the author investigated the emission of he which was found by CLEOc i 



n 



and 



intrigued great interests of theorists [69|] as well as experimentalists [70|]. The decay width 
of such transitions can be written as 

where |Ci p is a constant to be determined and it comes from the hadronization of gluons into 
pions, G is the phase space factor, f^fl'^^np ip overlapping integral over the concerned 
hadronic wave functions, their concrete forms were given in 671] as 

^ / RFir)rP^R*j^i{rydr J R*j^i{r'y^' Rj{r')r'^dr' 

K 



fl,Pl,PF 
J nj,lj,np,lf 



Mj - Eki 

where ni,np are the principal quantum numbers of initial and final states, Ii,If are the 
angular momenta of the initial and final states, / is the angular momentum of the color- 
octet qq in the intermediate state, Pj, Pp are the indices related to the multipole radiation, 
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for the El radiation Pi,Pf=1 and / = 1. Ri,Rf and Rki are the radial wave functions 
of the initial and final states, Mj is the mass of initial quarkonium and Eki is the energy 
eigenvalue of the intermediate hybrid state. 

The framework provides a more elegant way to deal with the long-distance QCD effects 
even though it only concerns transitions between states of heavy quarkonia. Moreover, 
between two gluon emissions the intermediate state is a hybrid state which definitely is a 
subject to draw interests of all high energy physicists. Just as aforementioned, the QCD the- 
ory predicts existence of exotic states, such as glueball, hybrid, tetraquark and pentaquark 
etc., at least does not exclude their existence, however, so far none of such exotic states 
have been experimentally identified yet, so that any direct or indirect information about the 
exotic states must be valuable. Since the subject is heavy quarkonia, the potential model 
is reasonable to describe their structures. Recently, many works are devoted to study the 
potential which can describe the hybrids which are composed of heavy quark, anti-quark 
and a gluon. In this case the quark and heavy quark reside in a color octet to keep the 
hybrid meson in a color singlet, so that the Coulomb potential between them is repulsive. 



The recent literature suggests a form by Swanson and Szczepaniak 71 1 

V{r) =br + -(l- 



r 



where the concerned parameters were given in 71[. Alternatively, Allen et al. proposed 



another potential form which includes a repulsive Coulomb piece as [7! 



K 



V{r) = - + J{bry + 2'Kb + Vo. 



where Vq is the zero-point energ y a nd other parameters were also depicted in 72 1 



When the works of Refs. 



66 



67| were done, there were not many data about the transi- 
tions available, so that the authors assumed that ?/'(4.03) is the ground state of hybrid \ccg) 
and accordingly obtained the concerned parameters in the potentials listed above. Recently, 
thanks to the great work of Belle, Barbar and CLEOc, much more data have been collected 
and they enable us to take an inverse strategy to study the problem. 

In our strategy, we treat the parameters in the hybrid potential as free parameters to 
be determined. Using all the available data on T{nS) T{mS) + vr + vr and tpinS) — 



ip{mS) + TT + TT as inputs, we apply the analysis with the form of defined in 73|| as 

j w^'' - w ry 

Jaw! 
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where i represents the i-th channel, W/'* is the theoretical prediction on the width of channel 
i, W^^^ is the corresponding experimentally measured value, AW^^^ is the experimental 
error. 

Carrying out all the procedures, we have obtained that the masses of ground states of 
the charmonium family and T family as 4.23 GeV (for charmonium) and 10.79 GeV (for 
bottonium). It is noted that they are not the physical states which are experimentally 
observed. It is quite understandable because the present knowledge may suggest that the 
gleballs and hybrids may not exist as real resonances with fixed masses and widths, but mix 
with hadrons with regular valence-quark compositions [74 1 and the physical states are the 
eigenstates of the mass matrices. 

However, from other aspects, these results are not accurate due to large uncertainties 
of the experimental data. Especially, when we reached the results, T{5S) has not been 
measured yet, and we did not include its transition to lower T members. 

Last year, the Belle Collaboration reported their measurements on T(5S') T(lS')7r"'"7r~ 
and T(55) T(25)7r+7r- with decays widths 0.59 ± 0.04(stat.) ± 0.09(syst.) MeV and 
0.85 ± 0.07(stat.) ± 0.16(syst.) MeV. These values are about two orders larger than the 
previously measured partial widths for dipion transitions between lower T resonances [tS]. 

Meng and Zhao suggested that the anomalous enhancement is due to the final state 
interaction [3]. Namely, because T(55') and T(45') are heavy and above the production 
threshold of BB, therefore they may first decay into a BB pair and then by a re-scattering, 
B — B would turn into T(mS') -|- tttt with m < 3. By this picture and fixing the concerned 
parameters within reasonable ranges, the enhancement may be understood. If it is true, 
one cannot further use the data of T{5S) — T{mS) + vrvr in our above calculations because 
the re-scattering contribution contaminates the whole picture and one is no longer enable 
to gain direct information about the hybrid intermediate state at all. 

In the calculation, we have found a strange phenomenon that for T{nS) — > T(mS')+7r+7r, 
the results are pretty stable, however for ilj{nS) — > tp{mS) +71 + 71 there exists a cancelation 
among large numbers with smaller numbers remaining. This is due to the closeness of the 
chrmonia masses and the hybrids, therefore the results on the charmonia transitions are not 
very reliable. 

On other aspect, such unstableness may be resulted in by a mistreatment of the charmonia 
transition. If the final state interaction is important as the authors of |3] suggested, for 
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bottonia transitions, it would also apply to the charmonia transition when masses of ip{nS) 
are above the production threshold of DD, taking into account such FSI effects, one may 
re-extract information about the direct transitions. As did for the potential model, one may 
expect that the parameters are universal for b and c cases and then we can reduce theoretical 
uncertainties in the calculation. 

Indeed, such information is very necessary for determining the model parameters and 
even judge the whole scenario of hybrids. Therefore we are expecting more data in the 
charmonia energy regions to be collected in Babar, Belle and even the LHCb as well as 
improvements of theory. 



VIII. THE X, Y AND Z RESONANCES 

The QCD theory and quark model have been proven to be very successful, however, it 
is by no means the end of the story. Both QCD and quark model have soft belly where 
many problems are not answered yet. Interesting, the two aspects in the quark model and 
QCD are connected to each other. In the QCD theory, thanks to the asymptotic freedom, 
the perturbation can be trustfully applied to evaluate any high energy processes, however, 
on other side, the low energy processes are governed by the non-perturbative QCD effects 
for which so far there lacks any reliable way to deal with. We have already briefly discussed 
this issue in previous section. Unfortunately, many real physics quantities are related to the 
low processes, such as the fragmentation in high energy collisions and hadronic transition 
form factors in hadron decays. On other aspect, the quark model demands that mesons are 
composed of a pair of quark and anti-quark, baryons consist of three quarks (or anti-quarks) 
and QCD interaction (i.e. strong interaction) binds all the constituents into hadrons. Both 
of the theories do not prohibit existence of exotic states, such as glueball, hybrids and multi- 
quark states (tetraquark, pentaquark etc.) or even favor their existence. The pentaquark 
was a hot topic for a while as several groups claimed that pentaquark containing an anti- 
strange quark or anti-charm quark were observed, then new data of most of the major 
labs gave completely negative results. Do they really exist or mix with regular hadrons as 



suggested in literature [77|? That is still an unsolved question. 

Recently the Babar, belle, CLEO and BES reported many newly observed resonances 
which are randomly named as X, Y and Z particles. In a review paper, Godfrey and Olsen 
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discussed this issue in some details 78|. There are many theoretical works devoting to the 
exciting field. 



A. Dlj{23l7), D,j{2m), D,j{2860) and Dsj{27l5) 



79 



80 



81 



whose spin-parity 



The discoveries of mesons D*j(2317) and Ds,j{2AQ0) 
structure are respectively 0"'' and l"*", have attracted great interests of both theorists and 
experimentalists of high energy physics, because they seem to be exotic. Bardeen et al. 



supposed that D*j{2317) and D,j(2460) are (0+, 1+) chiral partners of and D* [SSj i.e 



p-wave excited states of Dg and D* 



84| . By studying the mass spectra, Beveren and Rupp 



suggested that D*j{2317) and Dsj{2AQ0) are made of c and s |85|. With the QCD spectral 
sum rules, Narison calculated the masses of D*j{2317) and Dsj{24Q0) by assuming them 
as quark- ant iquark bound states and obtained results consistent with the experiment data 
within a wide error range 86|]. Very recently, considering the contribution of DK continuum 
in the QCD sum rules, Dai et al. obtained the mass of D*j{2317) which is consistent with 



experiments 



87|. Meanwhile, some authors suggested that L'*j(2317) and Dsj{24:60) may 



be of four-quark structure 



89 



90| 



Thus to clarify the mist of the structures of D*j(2317) and Z^sj(2460), serious theoretical 
works are needed. The studies of the productions and decays of D*j{2317) and Dsj{24Q0) 
are very interesting topics. Several groups have calculated the strong and radiative decay 
rates of D*j{2317) and Dsj{2460) in different theoretical approaches: the Light Cone QCD 
Sum Rules, Constituent Quark model. Vector Meson Dominant (VMD) ansatz , constit uent 



quark meson model, etc. 



91 



92 



93 



94 



95 



96 



97 



98 



991] . The authors of Ref . 



89 



l3 



im also 



calculated the rates based on the assumption that D*j(2317) and Dsj{24Q0) are in non-cs 
structures. Their predictions on the D*j(2317) and Z)sj(2460) decay rates are obviously a 
few orders larger than that obtained by assuming the two-quark structure. Recent Faessler 
et al. study the same subject assuming D*j(2317) as a DK molecule state using an effective 
Lagrangian approach 101 1. 

The semileptonic decay of is one of ideal platforms to study the productions of 
D*j{2317) and D,j(2460). Especially the Large Hadron Collider (LHC) will be run- 
ning in 2008, which can produce large amounts of Bg. ^Thus the measurements on 

Bs Dsj{2317,24Q0)W would be realistic. In Ref. 102l |. author calculated rate of 
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Bs 
of Ref . 



D 



2317 , 2460)Zi^ in terms of the QCD sum rules and HQET. Recently, authors 



103 



104| completed the calculations of Bs Dsj{2317, 2460)/z/ semileptonic decays 



in the QCD su m rules an d obtained large branching ratios. However, the results obtain ed by 



authors of Ref. 



103 



1041 ] are one order smaller than those given in Ref. 102| . In Ref. 105| , 



authors studied the same topic in terms of the Constituent Quark Meson (CQ M) rn odel. 

W e s tima ted by the authors of Ref. jl05 | and 



The branching ratios of 5, D,j(2317, 2460) 
that obtai ned i n terms of the QCD sum rules 

In Ref. 1061], authors use the heavy quark effective theory (HQET) and a non-relativistic 



103 



104j are of the same order of magnitude. 



model to evaluate the production rate of D*j{2317) in the decay of ip{4:A15), and find that 
it is sizable and may be observable at BES HI and CLEO, if it is a p-wave excited state of 
D,(1968). 

Because Ds.j(26 32) w as only obser ved by the SELEX collaboration 1071 ]. but not by 



Babar 



108| . Belle |l09l ] and FOCUS [11Q|, its existence is still in dispute, so we do not 



intend to include Dsj{2632) in this review. 

In the summer of 2006, the Babar collaboration observed a new cs state Dsj{2860) with 
a mass 2856.6 ± 1.5 ± 5.0 MeV and width T = (48 ± 7 ± 10) MeV. Babar observed it only in 
the D^K~^, D^Kg channel and found no evidence of D*^K ~^ an d D*~^Kg. Thus its quantum 



number should correspondingly is = 0^, 1 , 2+, 3 , ■ ■ ■ lll| . At the same time, the Belle 



collabora tion reported a broader cs state -Dsj(2715) with J = 1 in B'^ D D decay 
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113| . Its mass is 2715 ± llt^ MeV and width T = (115 ± 20: 



-36\ 
-32) 



MeV. 



According to the heavy quark effective field theory, heavy mesons form doublets. For 
example, we have one s-wave cs doublet (0~, 1~) = (-Ds(1965), -D*(2115)) and two p-wave 
doublets (0+,l+) = {D*j{2317), L>,j(2460)) and (l+,2+) = (D,i (2536), L>,2 (2573)) H- 
The two d-wave cs doublets (1~, 2~) and (2~, 3~) have not been observed yet. The possible 
quantum numbers of D,j(2860) may be 0+i2^Po), l~(l^/^i), l~i2^Si), 2+{2^P2), 2+{l^F2) 
and 3~{1^D3). The 2^P2 cs state is expected to lie around (2.95 ~ 3.0) GeV while the mass 
of the 1^F2 state will be much higher than 2.86 GeV. 



Dsj{28Q0) was prop osed as the first radial excitation of D*j{2317) ll4 . Ill5l ]. or as a 



J' 



3 cs state 



116| or as cs{2P) state 1171 ]. By the potential model, one can see that 



-Dsj(2715) sit s exa ctly at the position predicted by the quark model, as 2715 MeV if it is a 



2'^ 5*1 cs state 



is formed 119 ] 



1181 ]. The 1 state should lie around 2721 MeV if a (l"*", 1 ) cs chiral doublet 
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In Ref. 



120], authors investigated the strong decays of the excited cs states using the ^Pq 



model. After comparing the theoretical decay widths and decay patterns with the available 
experimental data, they tend to conclude: (1) Dsj(2715) is probably the l~{l^Di) cs state 
ahhough the l-{2^Si) assignment is not completely excluded; (2) Dsj{2860) seems unlikely 
to be the l~{2^Si) and candidate; (3) D,j(2860) as either a 0+(23Po) or S-ii'^D-s) 

cs state is consistent with the experimental data; (4) experimental search of Dsj(2860) in the 
channels D.^ y, D K*, D*K and D*?7 will be crucial to distinguish the above two possibilities. 



In Ref. 



121| . the strong decay of D wave cs meson to light pseudoscalar meson are 



studied in the framework of light cone QCD sum rule (LCQSR). 

Recently Dubynskiy and Voloshin proposed an interesting picture to explain the newly ob- 
served rich family of X, Y and Z. They suggested that the charmonium states, such as J 
ip{2S), rjci can be bound inside light hadronic matter, especially inside higher reson ance s 
made from light quarks and gluons and they named such states as hadro-charmonium 122 1. 
Definitely, this picture should undergo some more serious theoretical and experimental tests. 

Indeed, this is a wide world to be explored, which may help to testify the quark model 
and the low-energy QCD principles as well as all the working phenomenological models. 



IX. - MIXING 



This is an extremely interesting subject since a sizable mixing of D^ — D^ generally implies 
existence of new physics beyond the standard model. 



In the SM, mixing of partia 
occurs via the box- diagrams 



e an d anti-particle, such as — K", D^D^ and -Bj^-, — B[ 



mechanism 



123|. The calculation is standard based on the the GIM 



12J]. The contribution of the box diagram is proportional to mf/ 



w. 



here m,- 



is the mass of the exchanged quarks in the box and the CKM matrix elements 125l |. Thus 
for the — and Bg — Bs mixing, the exchanged quark is the top quark and the factor 
mf/rn^, > 1 becomes an enhancement, so that the mixing is large and such mixing has been 
reported to be observed long time ago and in the history it played an important role to hint 
that top quark might be heavier than the ly-boson. However, for — D^, the exchanged 
quark can only be b-quark which is much lighter than the top quark, so that the resultant 
mixi ng must be very small. The — mixing has indeed been measured by the Babar 
126| and Belle 1271 ] collaborations. Therefore it implies possible existence of new physics 
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beyond SM. There have been some theoretical suggestion about the mechanisms which may 
cause a siza ble — mixing, for exam ple, via a FCNC process in the non-universal Z' 
model 



128| or the unparticle model |129 |. 



The eigenstates of the mass matrix are 



130 



131 | 



\DH)=p\D')+q\D'), 

\DL)=p\D')-q\D°). 
with Ij^P + Igp = 1 and the corresponding eigenvalues are 

where M12 and are the off-diagonal matrix elements and obtained by calculating the box 
diagrams, it is noted that both of them are complex. In SM, they are small, thus a sizable 
non-zero mass and life difference of the two physical states \Ds) and |-Dl) must be caused 
by new physics, as aforementioned. 

Since the mass and life differences of the two eigenstates are not very large, it is hard to 
measure it as one did for the — K'^ system, we can investigate the evolution process of 
the -D° — system which in general is produced in S-decays or higher excited states of ijj 
family. One has 130 1 



\Dl{t))=g^{t)\D') + ^-g4t)m, 



g^{t)\D') + g^{t)m, 



and 



1 -jmt-^f 

9± = 2^ 



e 2 



A7j 



Am A7 , 



and Am = mn — ttll, A7 = '^h — 7l- The important parameters are 



Am 



X 



y 



A7 



7 7 
where 7 is the average lifetime of Dh and Dl. 

In analog to the i^-system, there exists the indirect CP violation which can be observed 
in the time evolution of the system. The direct CP violation will b e dis cussed later. 

The data of the Babar collaboration about the — mixing 126l | are x'^ = [—0.22 ± 
0.30(stat) ± 0.21{syst)] x 10"^ and y' = [9.7 ± 4.4(stat.) ± 3.1(syst.)] x 10"^ where x' = 
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X COS Sk-k + ysin^xTr; y' = — xsin^A'TT + y cos Sxn and Sk-k is the strong phase between the 
Cabibbo favored (CF) and doubly Cabibbo suppressed (DCS) amphtudes. x'is consi stent 



with zero and y' obviously deviates from zero. The data given by the Belle collaboration [127 1 
areycp = [1.31±0.32(stat.)±0.25(syst.)]% where y^p is defined as y^p = y cos (j)—^A^f sincf). 
When CP is conserved, Am = (f) = 0, the results are consistent with that obtained by the 
Babar collaboration. 

The unparticle model was first proposed by Georgi 1331 ]. Georgi argued that operators 



Obz made of BZ fields in the scale invariant sector may interact with operators Osm of 
dimension dsM made of Standard Model (SM) fields at some high energy scale by exchange 



particles of large masses, Mk, with the generic form OsmObz/M^. At another scale A, 



u 

the BZ sector induces dimensional transmutation, below that scale the BZ operator Obz 
matches onto unparticle operator Ou with dimension du and the unparticle interaction with 
SM particles at low energy has the form 

In the SM, the weak phase comes from the CP phase in the CKM matrix, since the SM 
contribution to M12 and can be neglected, the weak phase must be induced in the new 
physics. We would like to point out some salient features of the unparticle contribution to 
and due to an extra phase factor e"*'^'^'^ in the new scenario. We note that Mf^ 
can have both positive and negative signs depending on the value of du due to the factor 
cot(7r(i(/), therefore if information about the sign can be obtained from other theoretical 
considerations or experimental data, the dimension du would be restricted. There may be a 
sizeable contribution to at tree level which is not possible for usual mode where heavy 
particles are exchanged at tree level. For du equal to half integers, there is no contribution 
to M12, but there is to ri2. Another salient feature is that the unparticle contribution to 
the ratio Mi2/(ri2/2) is related to the unparticle dimension parameter du by 

<2 



ri'2/2 ~ 

If the unparticle contribution dominates meson and antimeson oscillation then the measure- 
ments of M12 and ri2 would provide a possible way to determine the dimension parameter 
du- In this case, we may gain valuable information about the dimension of unparticle du 
which so far cannot be directly obtained when mapping the operator Obz onto Ou and 
remains as a free and adjustable parameter in most of phenomenological research works. 
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Definitely, the contributions of new pliysics would be added to that of SM. The new 
physics effects also contribute to — and B^^^ — K^^^ mixing, however, in those cases 
the SM contribution obviously dominates over that from new physics and the effects of new 
physics would be smeared out as the measurements are not very accurate. By contrary, the 
SM contribution to D° - is negligible and all contribution comes from the new physics, 
thus measurements on it may provide us with an ideal place for gaining valuable information 
about new physics which is indeed the goal of all high-energy physicists. 

To make an accurate measurement, one needs a longer flight time before the D meson 
decays, thus the main labs to observe — would be the two -B-factories, and the LHCb 
which will be running soon will offer us another ideal spot to carry out such measurements. 
Even though for this observation, the BESlll does not have advantages for the kinematics 
because the linear momenta of the produced D mesons is small the relativistic time dilation 
does not apply, as the builder promised, the database will be very large and may greatly 
enhance the statistics, so that it m ay a" 



so be a possible lab to observe the mixing effects and 
explore new physics. Li and Yang 132[ | studied how to properly extract information about 
the mixing from data and make a efficient analysis. 

It is natural to ask if one can observe CP violation in /^-system. If Dh and are not 
CP eigenstates, there could be an indirect CP violation, however, since D decays faster and 
there are many channels available, one cannot determine the indirect CP violation as easy 

as in the K system by measuring rj^ and r/oo at different distances. 

The direct CP violation is defined as 

^ ... ^ r(D°(t) -> /) - r(D°(t) ^ /) 



130|. So 



which is a time-dependent measurable quantity. Some details are given in Ref. 
far, there is no report on the observation of CP violation at /^-system yet, even though the 
direct CP violation has been measured to be non-zero at i?-system. One may expect to 
make progress along the line in the future. 

X. CHARMED BARYONS 

Let us present our notations for the excited charmed baryons. Inside a charmed baryon 
there are one charm quark and two light quarks {u, d ot s). It belongs to either the symmetric 
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6f OT antisymmetric Sp flavor representation (see Fig. [T]). For the S-wave charmed baryons, 
the total color-flavor-spin wave function and color wave function must be symmetric and 
antisymmetric respectively. Hence the spin of the two light quarks is S=l for Gp or S=0 for 
3^. The angular momentum and parity of the S-wave charmed baryons are = or I"*" 
for 6f and = |^ for Sp. The S-wave charmed baryons are listed in Fig. [H where we 
use the star to denote |^ baryons and the prime to denote the = I"*" baryons in the 
representation. 



Ei*'++(™c) 




j:i''"(ddc) Sy"{dsc) n!."'"(ssc) S?(<isc) 



6 3 



FIG. 1: The SU(3) flavor multiplets of charmed baryons 

In Fig. [2] we present our notations and conventions for the P-wave charmed baryons. Ip 
is the orbital angular momentum between the two light quarks while Ix denotes the orbital 
angular momentum between the charm quark and the two light quark system. We use the 
prime to label the Hcj, baryons in the Gp representation and the tilde to discriminate the 
baryons with Ip = 1 from that with = 1. 

The notation for D-wave charmed baryons is more complicated (see Fig. [3]). Besides the 
prime, Ip and Ix deflned above, we use the hat and check to denote the charmed baryons 
with Ip = 2 and Ip = 1 respectively. For the baryons with Ip = 1 and /a = 1, we use the 
superscript L to denote the different total angular momentum in A^j^, S^j^ and '^cj,- 



The Babar and Belle collaborations observed several 



excited charmed baryons: 



Ae(2880,2940)+, 2^(2980,3077)+'° and ^^(2768)° [l34, ll35|, ll36|, ll37|, ll38|, ll39|. We col- 
lect the experimental information of these recently observed h adro ns in Table IIIII Their 
quantum numbers have not been determined except Ac(2880)+ 140 1. 
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'cO { 2 ) 
cl(j >J ) 
c-id > 5 ) 

) 

cl(5 .5 ) 

c2(i",r) 
ill I -I ) 



FIG. 2: The notations for P-wave charmed baryons. fsi^p) and fAi^p) denote the SU(3) flavor 
representation. Sg-^q^ is the total spin of the two light quarks. L denotes the total orbital angular 
momentum of charmed baryon system. 



State 


Mass and Width (MoV) 


Decay channels in experiments 


Otlicr information 


Ac (2880) + 


2881.5 ± 0.3, < 8 [1401 


Ac7r+ 77 ^ 


favors 1+ [1351. 
r(s;(2520),r ) ^ ^ [1351 


2881.9 ± 0.1 ± 0.5 , 5.8 ± 1.5 ± 1.1 [1341 




2881.2 ± 0.2+°-3, 5.5+° g ± 0.4 [1351 


E;"' + +(2520)ir+-" 


r(Sc(2455)7r±) 


Ac (2940) + 


2939. ± 1.3 ± 1.0, 17.5 ± 5.2 ± 5.9 [1341 






2937.9 ± l.olj- 5, 10 ± 4 ± 5 [1351 


Sc(2455)°' + +7r+'" 


Sc(2980)+ 


29G7.1 ± 1.9 ± 1.0, 23.6 ± 2.8 ± 1.3 [1361 


A+K"7r + 




2978.5 ± 2.1 ± 2.0, 43.5 ± 7.5 ± 7.0 [1371 


A+K"7r + 


Hc(2980)'' 


2977.1 ± 8.8 ± 3.5, 43.5 [1371 


A+Z^Ott- 




Sc(3077)+ 


3076.4 ± 0.7 ± 0.3, 6.2 ± 1.6 ± 0.5 [1361 


A+ K-TT + 




3076.7 ± 0.9 ± 0.5, 6.2 ± 1.2 ± 0.8 [1371 


A+K"7r + 


3^(3077)" 


3082.8 ± 1.8 ± 1.5, 5.2 ± 3.1 ± 1.8 [1371 


A+xn^r- 




Sc(3055)+ [1381 


3054.2 ± 1.2 ± 0.5, 17 ± G ± 11 


A+K-TT + 




Hc(3122)+ [1381 


3122.9 ± 1.3 ± 0.3, 4.4 ± 3.4 ± 1.7 


A+K-TT+ 




nc(2768)° 


2768.3 ± 3.0 [1391 




7P _ 3 + 

^2 



TABLE III: A summary of recently observed charmed baryons by Babar and Belle collaborations. 



(a) = 0. = 1 



/s(6): £ = 1®5„„ = 1< 



fA(3): L = l®5„„ =0 



Ji = 0: Scod ) 

,/, = !: Sci(l",r) 
J, = 2: Sc^dM") 

J, =l: Aci(l-,i-) 



(b) = 1, h = 



/,i(3): L = l®5„„ = 1< 



J, = 0: Acii(i ) 



/s(6); L = l®5„„ =0 



./, = !: ScidM ) 



In the past decades, there have been some research works on heavy baryons 
How ever these new observation inspired new investigations of these states 



141 



144 



142 



145 



431. 



146, 



1471]. In Ref. 145| . authors studied the A^940)"'" and its possible decay modes assuming 



Ac(2940)"^ to be a D*'^p molecular state [145|. Cheng et al. calculated the strong decays of 
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(a) Ip = 0,lx = 2 



/s(6): L = 2»S„„ = 1< 



/a(3): L = 2 » S,„3 = 



(b) lp = 2,l, = (} 



/s{6): L = 2CSS,,„ = 1< 



/a{3): L = 2 ® S„„ = 



(c) lp = l,h = l 



J, = 1: S,: 



J, =2: S,2(r,f^) 



=2: A,2(| + ,| + ) 



cll2 



1+ 3^ 
2 ' 2 



,3+ 5H 
-c2^2 ' 2 



V /5+ 7+ 

'cSlj '2 



- /3+ 5 + 



J, =2: S,2(i",r) 



.7, =2: A.2(f .i"-) 



--I /1+ 3^ 
-cl^2 '2 



; ,5+ 74 
'c3l2 '2 



=.c2(2 



/-4(3) 



/s(6) 



i = 0»5„,, = 1 



i = 1 ® S,,,^ = 



' J, = 0: AJ„(1-^) 

,7, = 1: AJ,(lM + ) 

-J, = 2: Ai2(iM^) 

^ = 1: A»,(lM + ) 

,J, = 1: A^,(lM + ) 

J, =2: A?2(iM^) 

J, =3: A?3(|M^) 



J, = 0: S°„(i + ) 
,/, = !: SjidM^ 

J, =2: s?2(r,r 



"Jo(j ) 



"e2(,2 



"elV2 



^•2 ,-3+ 5+\ 
"e2(,2 '2 ^ 



S2 ^5+ 74 
"e3(,2 '2 



"c()(.2 J 

3 + 
"cl^2 '2 

^/2 / 3 
"<:2t2 



FIG. 3: The notations for the D-wave charmed baryons. 



newly obser ved c harmed mesons in the framework of heavy hadron chiral perturbation theory 
(HHChPT) 1461]. I n order to understand their structures using the present experimental 
information, in Ref . 148l . Il49| , the strong decay pattern of the excited charmed baryons are 
studied systematically in the framework of the ^Pq strong decay model. After comparing the 
theoretical results with the available experimental data, their favorable quantum numbers 
and assignments are obtained in the quark model. 

It is an interesting field. There was a puzzle that the lifetime of is close to that of 



B-- 



r 



(5" 



: (1.530 ± 0.009) X 10-12 ^(^±) ^ ± g.Oll) x lO'^^ whereas the 

lifetime of is almost double of that of {t{D^) = (410.1 ± 1.5) x 10-^^ s, r{D^) = 
(1040 ± 7) X 10"^^ s). By the heavy quark effective theory, the total width of a meson is 
mainly determined by the total width of the heavy quark constituent in the meson which 
is proportional to rriQ x Vckm where Vckm is the corresponding CKM matrix elements. 
Obviously, the weak decays of b-quark is Cabibbo suppressed while for charm quark, it is 
Cabibbo favo red. The obvious difference between B and D lifetimes is excellently explained 
by Bigi et al. 150l ]. Another puzzle is the lifetime difference between A;, and B meson, which 
used to be as large as 0.29, howe ver the recent measurement alleviates the disc repancy as 
T{Ah)/T{B^) = 1.041 ±0.057 1511 ]. and is close to the theoretical evaluation 1521 ]. However, 



in the theoretical calculation one needs to evaluate the hadronic matrix element in analog 
to the discussion given in previous sections about the meson case, and the estimate still has 
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a large uncertainty, so we do not think that the problem is completely solved. 

More intriguing issue is that the lifetime of r(Ac)(~ (200 ± 6) x 10~^^ s) is much smaller 
than that of and D^, and it implies that the light quark cannot completely be treated 
as spectators and the binding effects may also be important. In the heavy quark language, 
the ^ {n > 1) corrections may play important roles when the total width is evaluated. The 
question is that for exclusive decays, can one use the pQCD or similar method where the 
l/mg exp ansio n is adopted? It seems that one can apply the pQCD method to analyze the 
Af, decays 153l |. but so far, there is not work devoting to pQCD application in calculating Aj, 
decays yet. The reason is that the energy scale of is relatively low and the exchanged 
gluons in the processes are not hard enough to enable us to consider only the leading, at 
most the NLO contributions. How to involve more higher terms in the expansion may be a 
very challenging task for theorists. 

A great progress in the fiel d is that the double- charmed baryon Sec was observed by SE- 



LEX 



154 



155 



156 



157, 



158l | at the Fermi Lab. Since it contains two charm quarks which 



would decay via weak interaction independently, one may expect to gain some information 
abou t the binding effects. In analog with the method given by Bigi et al. 1501 ] . we evalu- 



ated 



159| the lifetimes of 



' — 'cc^ ' — 'cc 



and which contain different light flavors, to order 
of l/m^. The hadronic matrix elements were evaluated in terms of the simplest model, i.e. 
the harmonic oscillator model 36|]. The theoretical predictions can be tested in the future 
experiments especially the LHCb. Indeed, we lay our hope on the LHCb data which may 
produce a large database on such doubly charmed baryons, and maybe also baryons con- 
taining two b-quarks H^b, or b and c quarks S^c etc. By studying them, we will gain valuable 
information on the non-perturbative QCD effects which bind the heavy quarks together with 
a light flavor. Moreover, one may expect to observe baryons which are compose d of three 
heavy quarks {b and/or c). There has been some theoretical works in this field |l60 |. but 
phenomenological studies on their production rates and decay modes require more serious 
works in order to investigate their structures and governing mechanisms in some details. 
Some details of baryons containing two heavy quarks have been reviewed by Kiselev and 



Likholded in an enlightening paper 



161 |. 



We have also proposed a special decay mode of H^c to investigate new p hysic s. Namely, 



we suggest to measure a direct decay of S^c into non-charm final states 



162l |. which is 



definitely a signal for new physics beyond SM. Concretely, in the SM, the direct decay 
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of Sec into non-charm final states are very restricted and tfie sequential decays, such as 
Sec ^Ac + -D^A + A' + X would produce final states which involve at least more than 
two hadrons, so that can be easily distinguished from the expected decay modes with two 
non-charmed hadrons. We have tried to employ two models which recently become more 
popular to the theorists, the Z' and unparticle models which can induce the flavor-changing 
neutral current at tree level, so can realize the reaction cc — > gg' where q and g' are light 
flavors. However, our numerical results show that the two models cannot cause sizable 
fractions for practical measurement, even though one can have a great luminosity at LHCb. 
Our conclusion is that if such decay modes are observed, one can claim that as a definite 
evidence of new physics, but it is neither the Z' model, nor the unparticle model. This 
would motivate us to explore other new physics models beyond the SM. 

The charmed pentaquark was a hot topic and the QCD theory does not exclude pen- 
taquark, the question is do we need to take it more seriou sly a nd where should we search for 
them and how to identify them from the regular baryons 163l ]. As we discussed in previous 
sections, many theorists favor tetraquark as a plausible explanation of some newly observed 
X, Y and Z resonances, therefore why should we firmly reject existence of pentaquark? A 
full scan of the pentaquark is not intended in this review work, but we would like to analyze 
some aspects concerning the charm physics. The first proposal on existence of pentaquark 



is based on the group theory analysis 



164| . that the pentquark is composed of four u and 



d quarks and one anti-strange quark with both the strangeness baryon number being 
Then a claim of finding charmed pentaquark was made where a c replaced s. However, later 
most major labs in the world reported negative results in search for pentaquark and the 
society of high energy physics tends to deny discovery of pentaquark. We also made effort 
to look for trace of pentaquark via some processes, such as photo-production and radiative 



decay of pentaquark as well as the hadronic decays involving pentaquarks [165j, and the 
conclusion is still that more accurate measurements are necessary. On other aspect, we 
seriously consider that the pentaquark may mix with the other baryonic states 166 |. 

It is interesting to investigate the quark-structure of pentaquark which is favored by 
QCD. Besides assuming that all the four quarks and an anti-quark all mix together in a big 
hadronic bag, the mor e favorable structures are proposed by several authors, for example. 



Jaffe and Wi 



cze 



167l | suggested the diquark-diquark-anti-quark structure whereas Karliner 



and Lipkin 168l | favored the diqaurk-triquark configuration. The key point is why the 
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pentaquark has not been observed in any colliders but was seen in fixed target experiments. 
Lipkin tried to explain the situation by considering the initial quark configuration in the 
beam particle or target. This should draw serious attentions of experimentalists and theorists 
and design new experiments to search for evidence of pentaquarks. Since in the pp collisions 
at LHC, the rich quark contents and high luminosity may be a good source to produce 
pentaquarks, and until then we can draw a definite conclusion if the pentaquark can exist as a 
real particle. In this line the charmed pentaquark may be more favorable because it contains 
a heavy fiavor. The most interesting subject would be that if the pentaquark is strictly 
prohibited, there must be a symmetry to restrict it, it would require a new understanding 
on the mechanism beside the general theory of QCD. 

XI. DIQUARK 

It is an extremely interesting topic, not only for charmed baryons, but since it is widely 
applied in studying physical processes where baryons are involved, diquark is worth careful 
investigation. By the SU(3) theory, two quarks can reside in a color- anti-trip let 3 to consti- 
tute a loose bound state and then it combines with the rest quark to make the baryon of 
color singlet. The QCD induced potential is proportional to the Casimir factor {iIj\X°'X°'\iIj) 
where A" is the SU(3) generator and is the state, thus in a color-anti-triplet (or a triplet 
for two anti-quarks) the two quarks attract each other because the Casimir factor is nega- 
tive. It seems that the bound state is plausible, but for light quarks, their linear momenta in 
the diquark are rather large and may spread out in space, thus cannot make a real physical 
subject in common sense. Therefore whether the diquark structure is only a mathematical 
description for a baryon or it can be a physical subject is still in dispute. That is how to 
understand the subject "diquark" is an intriguing topic in theoretical physics. 

However, if the baryon contains a heavy quark, it becomes simpler as usually assumed, 
the heavy quark may sit near the geometric center of the hadron (almost at the center) and 
one can further postulate that the linear momentum of the heavy quark is zero. Then the 
two light quarks would compose a diquark even though its spatial spreading is still large. 
Therefore when we apply the concept of diquark in this case, we indeed do not treat it as a 
point-like particle, but a loosely bound state with a common momentum and interact with 
the heavy quark as a whole object, moreover, when some reaction takes place, it may also 
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perform as a whole object, but to compensate the effects of its inner motion, one needs 
to introduce a form factor(s) whose phenomenological form was given in Ref.anselminol as 
F(Q'^ ) = q?^°q9 , where Qq is a parameter and is determined as Ql ~ 3.2 GeV^ by fitting 
data 169| . The form is obviously understandable. The form factors should be normalized 
to unity as 0, i.e. as one looks at the diquark from a far distance, the form factor 

becomes a unity, whereas a.s Q"^ oo, the inspector then penetrates into the diquark, so that 
he would see the individual quarks instead of the whole and the diquark picture no longer 
holds and mathematically it is required to approach zero as Q"^ ^ oo. The form factors 
should be calculated in terms of a more fundamental way, i.e. based on quantum field theory. 
However, the non-perturbative QCD effects can by no means be treated completely based 
on an underlying theory so far, instead we need to invoke phenomenological models. The 
Bethe-Salpeter equation obviously is a reasonable approach where the kernel depends on 
concrete models. The form factor obtained in terms of the BSE generally coincides with 



the picture described above. Guo, Thomas and Williams |l7Cll | studied the l/mg corrections 



for the B-S equations for Ag, Qq in the diquark picture. Then in the fram ework of the 



B-S equation, we calculate the form factors for various types of diquarks 17l| and the form 



'acto rs obtained based on the B-S theory are qualitatively consistent with that given in Ref. 



169| . Appl ying the diquark picture, we calculate several processes where heavy baryons 



participate |172| . The results are somehow competent to be compared with data, and it 
indicates that the diquark picture is vigorous and robust, even though still needs further 
studies. 

It may be worth mentioning a special situation, namely the bary ons containing two 



heavy quarks should fit in the diquark picture well. Falk et al. |l73l ] considered the two 
charm quarks in a doubly heavy baryon constitute a "perturbatively bound diquark" whose 
wavefunction at origin follows 

|i?(,,)(0)p^^|i?V^(0)|2^(0.41 GeV)^ 

and the factor 1/8 is due to the difference of the color structures of color anti-triplet and 



smg. 
try 



.et_. Georgi and Wise proposed a very special symmetry called as the superflavor symme- 
Simulating the supersymmetry where the SM fermions have their supersymmetric 



scalar partners and the SM bosons also have their SUSY fermionic partners, in the superfla- 
vor symmetry, the heavy quark is written in a doublet with its superflavor partner (scalar 
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or vector) of color triplet as 
and 

where is a heavy quark and obeys = h^, Xv and A'^ are heavy color triplet scalar 
and vector respectively, with constraint = 0. With this symmetry, we can estimate the 
production rates of heavy baryons which include two heavy quarks, i.e. a heavy diquark. In 
this scenario, the heavy scalar or vector diquarks can be treated as the superflavor symmetric 
partners of heavy quark and then one can compare the baryon processes with corresponding 
meson cases where two heavy mesons are produced. For both cases of baryon and meson, 
the light quark would play the same role. In the meson cases, at the leading order, there 
is only one Isgur-Wise function to manifest the non-perturbative QCD effects, thus by the 



superflavor symmetry, we would also apply the same function [28|] to calculate the production 
rates of Xcc, X^c, ^feb, etc. in e~^e~ collisions. Definitely the same procedure can be applied 
to calculate the heavy baryon production rates at LHC, even though the situation is a bit 
more complicated. 

XII. PRODUCTION OF CHARM 

Production of charm, i.e. production of charmed mesons, baryons and charmonia, is a 
large subject, and we are not going to cover this field in this short review, but indicate some 
tricky and challenging issues. 



The subject of production of J/ip, ip' and T was reviewed by Lansberg 175| in some 
details. In this incomplete review, we do not intend to cover the whole subject, but indicate 
that this field is still tricky and intriguing by two examples which are widely discussed in 
the society. 

As a long known question, the inclusive J/ip production in e^e~ J/tpcc raises a 
challenge to our theory. At LO, the QCD prediction on the cross section is smaller than 
the observed value by a factor of 5. This discrepancy may come from the way efunctions of 



the produced J/ip 174| or the NLO corrections. In an enlightening work 1771 ] . the authors 
suggested the NLO may enhance the cross section by a factor about 1.8. This may greatly 
change the situation. 
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However another serious question is raised that if the NLO can change this cross section 
so much, what mechanism or symmetry determine this unusual enhancement. Generally the 
loop suppression is proportional to a^/ vr and at this energy scale as is smaller than 0.2. If the 
number of diagrams (at NLO, the number of diagrams is rather large) is the reason as long 
as the interference among the diagrams are mostly constructive, one could ask how about 
the contribution of NNLO? Definitely, even though the diagrams at NLO are constructive, 
it does not suggest that they are constructive at NNLO at all. If the NNLO contribution 
decreases, we can happily declare that the theory works well, however if it continues to 
increase or the correction at NNLO is even larger than that of NLO, the unitarity would 
confront a serious challenge and one needs to take the hint more seriously. Indeed, to give a 
firm answer to the question, one must calculate the NNLO corrections. However, there are 
over 1000 two-loop diagrams and a complete calculation is extremely difficult project. Since 
it is necessary, some serious work must be done and it is indeed worthwhile. 



The problem about the production of J/ij) was investigated by several groups |l76l |. The 
authors studied the hadronprocution of J/V^, V^', T in associated with heavy quark pair 
which was the subject we listed above and investigated further, and a possible solu tion 
of the J 111) pro duction puzzle. Their result was confirmed by Gong and Wang 1781 ]. In 



the paper [178|, the authors further studied the pt distribution of the J/i}) polarization in 
pp ^ J /ip + X process. There is also serious discrepancy between theoretical prediction and 
data. They calculated the NLO pt distribution of J/ip polarization and found an amazing 
change. Namely, they find that the polarization status of J / ip at large pt changes from the 
transverse-polarization dominance at LO into the longitudinal polarization dominance at 
NLO. This correction indeed partly makes up the gap between theoretical result and data. 
In their first work, they only considered the color-singlet contribution and got improvement, 
so there was a hope that when the color-octet contribution was involved the situation would 
be improved further and the data could be explained eventually. But unfortunately, when the 
color-octet contribution is taken into account, the new correction is not so drastic and does 
not change the whole picture muc h, th erefore the large gap between theoretical prediction 



and data still exists. The authors 



178l | claimed that to explain the data, new mechanisms 



or new physics or new understanding of the J/ip structure etc. are needed. 

The two examples do not intend to cover the large-scale field of production of charm, 
but indicate that there are still many problems which are not fully understood yet, a lot of 
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theoretical efforts are required. This is also associated with the puzzle we discussed in this 
review that why pentaquark has not been observed and so on. 



XIII. DISCUSSION 

This is by no means a definite conclusion because in this review, we raise many unsolved 
problems and list some references where many authors devote great efforts in this fascinating 
field. The theoretical difficulty is obvious that charm quark is heavy, but not sufficiently 
heavy, so non-perturbative QCD effects still play important roles, and moreover, the rela- 
tivistic effects of quark quark in hadrons are also not negligible. Therefore the heavy quark 
effective theory (HQET) may apply to this field, but needs to consider higher orders in the 
l/rric expansions and it brings unexpected uncertainties. The good example is the lifetime 
difference between and D mesons, compared to the 5-cases. On other aspect, careful 
studies in the charm field may greatly enrich our understanding of the basic principles, 
concretely, how to properly deal with the non-perturb ativ e QC D effects is one of the most 



important issues. Moreover, the work by the authors |177l . Il78l ] indicates that higher order 
QCD corrections may be very important, and even play a crucial role. This indeed further 
intrigues the field of charm physics. 

The main task in this field seems not closely related to searching for new physics beyond 
SM, as generally considered because of its energy scale, but there indeed is possibility to 
investigate new physics. For example, the — mixing as aforementioned, and the 
lepton flavor violation processe s be sides the neutrino oscillation experiments. Namely, one 
can examine the reactions like 179[ | 



e+e" ^i)^ e-(+V^^"\ 

for which the charm-tau energy region may be the ideal observation place. On the theoretical 
aspect, in the SM, this FCNC can occur at one- loop order where neutrino masses need to 
be large enough to produce observational effects. However, it contradicts to all the data of 
the solar neutrino, atmospheric neutrino experiments and one can claim that observation of 
such reaction is a clear signal of new physics beyond SM. There are so many new physics 
models which can produce this reaction, can the reaction, if it indeed occurs, distinguish 
among them? Generally it cannot. But one can at least tell which model is a possible one 
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and maybe, such a reaction would help to get rid of a few proposed models. 

There are gen eral discussions on possibilities to measure the production rates from the 
;heo retical angle 18Cll | and set bounds on the model parameters. Very recently, we suggest 



181| that the unparticle model may induce such FCNC. By adopting the reasonable region 



of the concerned parameters, our calculation shows that one can expect a large production 
rate near the resonance peak of J/ip. Carefully analyzing possible background, one has 
confidence to extract clear signal for the luminosity of PEPC II and detection accuracy of 
BESIII. If it is observed, the unparticle model at least is one of favorable candidates of 
new physics, on other side, if it is not observed, the unparticle model would be somehow 
disfavored. 

There are too many questions to be answered in this energy region, as we discussed 
above ths,^=^.al ve.y enUgMening .view pape. which =hed l.gM. o„ th. excit.g 



field |7l.ll82l.ll83l.ll84|. 



Indeed, besides the major line, there are some other subjects which should be researched 
in this field, for example , it i s proposed to test the Bell inequalit y in h igh energy processes, 
Tjc AA — s> p-K^ + pTi^ I85I and decays of charmonia into kaon 186 |. As a matter of fact, 
the scheme proposed in 185l | was carried out in terms of the DM2 data, however, because 
the database was too small, the poor statistics made the work meaningless. Today, the 
situation will be different, the BESIII will provide a tremendously large database of J/^/', so 
that one may expect to make accurate measurements to testify the Bell inequality. 

We believe that we are standing at a prosperous epoch for high energy physics. The 
two B-factories continue to produce more data where one may also have a chance to study 
charm-related physics, such as CP violation and higher excited states of charmonia and 
charmed mesons and baryons. Moreover, the LHC will be running soon and not only 
a lot of information towards underlying physics such as the Higgs mechanism and new 
physics beyond SM will be collected, but also some details about bottom and charm 
physics can be achieved, especially we may lay great hope on the LHCb. Even the 
ALICE can tell us some interesting subjects about charm physics, for example, if the 
J I'll) suppression in quark-gluon-plasma (QGP) really exists, can the higher temperature 
and pressure change the potential between quarks, can we expect an observable phase 
transition, like deconfinement and chiral symmetry restoration etc. More than anything 
else which can enrich our knowledge on charm physics, the BEPC II and BESIII will 
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be operating very soon, and a great database on J/ip and other members of the iIj 
family, as well as the D-mesons, baryons, even tensors will be available, so that we are 
sure, may new challenges will be waiting for us. So we feel much encouraged than any 
past time, because so many questions need to be solved and more data will help to do the job. 
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